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Abstract 

A detailed analysis of available in situ and remotely sensed N 2 O and CH 4 data measured 
in the 1999-2000 winter Arctic vortex has been performed in order to quantify the temporal 
evolution of vortex descent. Differences in potential temperature (0) among balloon and aircraft 
vertical profiles (an average of 19-23 K on a given N 2 0 or CH 4 isopleth) indicated significant 
vortex inhomogeneity in late fall as compared with late winter profiles. A composite fall vortex 
profile was constructed for November 26, 1999, whose error bars encompassed the observed 
variability. High-latitude, extravortex profiles measured in different years and seasons revealed 
substantial variability in N 2 0 and CH 4 on 9 surfaces, but all were clearly distinguishable from the 
first vortex profiles measured in late fall 1 999. From these extravortex- vortex differences, we 
inferred descent prior to November 26: 397±15 K (la) at 30 ppbv N 2 0 and 640 ppbv CH 4 , and 
28±13 K above 200 ppbv N 2 0 and 1280 ppbv CH 4 . Changes in 0 were determined on five N 2 0 
and CH 4 isopleths from November 26 through March 12, and descent rates were calculated on 
each N 2 0 isopleth for several time intervals. The maximum descent rates were seen between 
November 26 and January 27: 0.82±0.20 K/day averaged over 50-250 ppbv N 2 0. By late winter 
(February 26-March 12), the average rate had decreased to 0.10±0.25 K/day. Descent rates also 
decreased with increasing N 2 0; the winter average (November 26-March 5) descent rate varied 
from 0.75±0.10 K/day at 50 ppbv to 0.40±0.1 1 K/day at 250 ppbv. Comparison of these results 
with observations and models of descent in prior years showed very good overall agreement. 

Two models of the 1999-2000 vortex descent, SLIMCAT and REPROBUS, despite 0 offsets 
with respect to observed profiles of up to 20 K on most tracer isopleths, produced descent rates 
that agreed very favorably with the inferred rates from observation. 
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1. Introduction 

Each fall, decreasing polar sunlight triggers a series of events resulting in increased 
circumpolar winds, the formation of an isolated region of stratospheric air known as the polar 
vortex, and the subsidence or diabatic descent of the vortex air column by several km resulting 
from radiative cooling [Schoeberl and Hartmann, 1991; McIntyre, 1992]. Quantifying vortex 
descent is necessary if one is to understand chemical processes, such as ozone loss, occurring 
within the vortex, as it is otherwise difficult to differentiate between chemical and dynamical 
changes in ozone. 

Quantitative studies of vortex descent have been carried out in the past. Tracers such as 
N 2 O, CH 4 , HF, CO, HCN, and CCI 2 F 2 have been used to tag similar air masses when comparing 
measurements made at different times and locations, and from these differences, descent has been 
inferred in the Arctic [Loewenstein et al., 1990; Mankin et al., 1990; Schmidt et al., 1991, 
Schoeberl et al, 1992; Toon et al, 1992; Bauer et al, 1994; Emmons et al, 1994; Traub et al, 
1995; Abrams et al, 1996a; Muller et al, 1996, 1997; Hartmann et al, 1997; Randel et al, 1998; 
Manney et al, 1999] and Antarctic [Jaramillo et al, 1989; Loewenstein et al, 1989; Toon et al, 
1989; Schoeberl et al, 1992, 1995; Russell et al, 1993a; Tuck et al, 1993; Crewel! et al, 1995, 
Abrams et al, 1996b; Randel et al, 1998; Manney et al, 1999; Allen et al, 2000; Kawamoto and 
Shiotani, 2000], A number of models have also been used to estimate descent [e.g., Schoeberl et 
al, 1989, 1990, 1992 \ Fisher et al, 1993; Manney et al, 1994, 1995, 1999; Rosenfield et al, 

1994; Strahan et al., 1994, 1996; Bacmeister et al, 1995; Eluszkiewicz et al., 1995; Lucic et al, 
1999]. Models have many advantages, including the ability to simulate the stratosphere in areas 
where measurements are unavailable. However, it is essential that such models be compared with 
observational results to validate their performance. 
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Satellite observations are extremely useful for this goal, as satellites offer near-continuous 
measurement and coverage over significant portions of the Earth. However, they also suffer from 
some major drawbacks, such as low horizontal (> 100 km) and vertical (~2 km) resolution, and 
for almost all solar occultati on-based instruments, limited coverage at high latitudes during the 
winter months. While the horizontal resolution is not particularly important for characterizing 
descent, which affects the vortex as a whole, the vertical resolution can be comparable to the 
amount of descent for short (< ~1 month) time periods. The inability to view polar latitudes 
during winter also hampers observations inside the vortex, which is typically located poleward of 
60 N or S. In the Arctic, stratospheric conditions sometimes displace the vortex away from the 
pole during portions of the winter, allowing plentiful observations inside the vortex [e.g., Muller 
et al., 1997; Manney et al., 1999]. However, these conditions are not necessarily frequent in each 
winter, as was the case during much of 1999-2000. 

Field-based observations (consisting of aircraft- and balloon-based in situ and remote 
sensing measurements, as well as ground-based remote observations) offer the advantages over 
satellite observations of high spatial resolution, especially for in situ measurements, and greater 
vortex penetration in winter months. However, they also suffer from several limitations, namely 
a lack of temporal coverage (field campaigns are expensive, and weather limits both balloon and 
aircraft frequency), spatial coverage (there are a limited number of polar locations suitable for 
airborne research), and altitude (the ER-2 aircraft is limited to ~20 km, and balloons are limited to 
~30 km). Thus, there is no instrument/platform combination that provides ideal vortex coverage. 

In order to study descent within the Arctic vortex during the 1999-2000 winter season, 
we have combined tracer observations from ten different instruments deployed as part of the 
Stratospheric Aerosol and Gas Experiment III (SAGE III) Ozone Loss and Validation Experiment 
(SOLVE) and the Third European Stratospheric Experiment on Ozone 2000 (THESE02000) 
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campaigns. These measurements represent a variety of techniques and covered large portions of 
the vortex from late November 1999 through mid-March 2000. In addition, measurements from 
three other instruments deployed in prior years were examined to help determine a vortex starting 
profile. Our inferred descent rates were compared with those from previous studies of Arctic 
descent. Lastly, our descent rates were compared with results from two stratospheric chemical 
transport models that simulated vortex descent for 1999-2000. 


2. Data and Methods 

In this paper, we focus on the stratospheric tracers N2O and CH4, which were measured 
by a large number of instruments during the SOLVE and THESE02000 campaigns. These tracers 
have long photochemical lifetimes (years to decades) in the middle and lower stratosphere and 
thus are ideal for studies of polar descent within the stratosphere, where transport timescales are 
on the order of months to years [Brasseur and Solomon, 1995]. 

Potential temperature (0) was chosen as the vertical coordinate because, unlike altitude or 
pressure, it is conserved under adiabatic displacements. In the absence of mixing, changes in 0 are 
due exclusively to diabatic heating, which is the effect we seek to quantify. To calculate 0, 
measurements of temperature and pressure are needed. In general, both these quantities were 
measured locally on the same platform as the tracer instrument; however, for some data products, 
assimilated pressure and temperature fields were used to generate 0. 

2.1. Instruments 

The instruments used in this study are listed in Table 1, along with platform, launch 
location, utilized data (N2O, CH4, 0), and flight dates. 
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The in situ balloon instruments will be discussed first. The Lightweight Airborne 
Chromatograph Experiment (LACE) determined the volume mixing ratios (VMRs) of eight gases 
including N 2 0 by gas chromatography (GC) with a time resolution of 70 s [Moore et al., 2001]. 
Pressure and temperature associated with the LACE flights were measured by the Jet Propulsion 
Laboratory (JPL) Ozone instrument [, Salawitch etal., 2001] on November 19, 1999, and by the 
National Oceanic and Atmospheric Administration (NO A A) Frost-Point Hygrometer [ Vomel et 
al., 1995] on March 5, 2000. The J. W. Goethe-Universitat cryosampler (BONBON) [Schmidt et 
al., 1987], which was located at the Forschungzentrum Jiilich prior to 1996, obtained up to 16 air 
samples per flight; VMRs were determined via GC in the laboratory. Pressure and temperature 
associated with the BONBON flights were measured by sensors located on the balloon payload, 
and operated by the Centre National d'Etudes Spatiales (CNES) and J. W. Goethe- 
Universitat/Forschungzentrum Jiilich teams, respectively. 

Data from the four in situ ER-2 instruments that measured N 2 0 were combined into a 
unified N 2 0 (uN 2 0) data product, following a procedure described in Hurst et al. [2001]. This 
procedure produced a self-consistent, high time-resolution (3 s) N 2 0 dataset for each SOLVE 
flight, using an objective method to evaluate and reduce bias among the instruments, followed by 
weighted averaging. The Argus spectrometer [Jost and Loewenstein, 1999] measured the VMRs 
of N 2 0 and CH 4 using second harmonic detection of an infrared tunable diode laser (TDL) in a 
multipass absorption cell, with an average time resolution of 3.2 s. The Aircraft Laser Infrared 
Absorption Spectrometer (ALIAS) [Webster et al., 1994] is a TDL spectrometer that measured 
the VMRs of four gases, including N 2 0 and CH 4 , with an average time resolution of 1.4 s. The 
Airborne Chromatograph for Atmospheric Trace Species IV (AC ATS IV) [Elkins et al., 1996; 
Romashkin et al., 2001] determined the VMRs of eleven gases, including N 2 0 and CH 4 , by GC 
with time resolution of 70 s. The Whole Air Sampler (WAS) [Schauffler et al., 1999] filled up to 
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32 canisters per flight with ambient air, which were analyzed in the laboratory for ~27 trace gases 
including N 2 O and CH 4 via GC. Pressure and temperature were measured on the ER-2 by the 
Meteorological Measurement System (MMS) instrument [ Scott et al., 1990]. 

Next we will discuss the remote sensing balloon instruments. The MkIV instrument is a 
Fourier Transform Infrared (FTIR) spectrometer, which measured solar occultation spectra over 
the wavelength range 1.77-15.4 pm. A large number of atmospheric constituents were measured, 
including N 2 O and CH 4 . Vertical resolution was 2 km, with profiles reported on a 1 -km vertical 
grid [Toon, 1991]. Pressure and temperature profiles, respectively, were derived from spectral 
fits to temperature-insensitive and temperature-sensitive CO 2 lines, assuming the CO 2 profile 
measured by the Harvard CO 2 instrument on board the Observations from the Middle 
Stratosphere (OMS) in situ gondola launched on November 19, 1999. 

The Submillimeterwave Limb Sounder (SLS) instrument recorded emission spectra of 
several gases; N 2 0 was measured near 600 GHz. The data presented here were flight averages of 
one or more limb scan emission spectra. Profiles were obtained from limb scan spectra using an 
iterative least-squares fitting procedure, with a nominal vertical resolution of 2.5 km. 

Temperature and pressure profiles were obtained from a balloon sonde launched on the day 
before the flight [Stachnik et ai, 1992]. 

On board the DC -8 aircraft was the Airborne Submillimeter Radiometer (ASUR), which 
retrieved complete vertical profiles of several gases, including N 2 O, every -2000-6000 s from 
pressure-broadened emission lines at 652.833 GHz. The vertical resolution was 5-10 km, with 
profiles calculated on a 2-km vertical grid. Flight track data containing temperature and pressure 
for calculating 0 were provided by the NASA Goddard Data Assimilation Office (DAO) [von 
Konig et al., 2000], 

Data from two spacebome instruments are used in this paper. The Halogen Occultation 
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Experiment (HALOE) instrument is located on the Upper Atmosphere Research Satellite 
(UARS), and uses solar occultation to measure vertical profiles of several gases, including CH4. 
The data used here were level 2, version 1 9. The instrument vertical resolution was 1 .6 km at the 
Earth’s limb, and the altitude range extended from 15 km to more than 60 km. Latitudinal 
coverage was from 80°S to 80°N over the course of one year. Because of the HALOE orbit and 
the Earth-Sun geometry during the Arctic winter months, the maximum northern coverage was 
limited to ~-50“N from November 1999 through January 2000, and up to ~65°N in October 1999 
and February-March 2000. As a result, data throughout most of the winter were limited to 
extravortex regions. Temperature and pressure were determined by spectral fitting of the 2.8 |xm 
CO2 band [ Russell et al, 1993b] above 35 km (10 mbar); below this altitude, they were obtained 
via interpolation from the National Meteorological Center (NMC) assimilated data product. 

The Atmospheric Trace Molecule Spectroscopy (ATMOS) instrument is an FTIR 
spectrometer similar to MkIV, which has flown on several Space Shuttle missions. ATMOS 
measured the abundances of ~30 gases, including N 2 0 and CH4, from 12-80 km altitude with a 
vertical resolution of 2-3 km. Version 3 data were used here. Northern latitude coverage ranged up 
to 69°N during the second Atmospheric Laboratory for Applications and Science (ATLAS-2) 
mission (April 8-16, 1993), and up to 49°N during the ATLAS-3 mission (November 3-12, 1994) 
[Gunson et al, 1996], Pressure and temperature were determined by spectral fitting of the C0 2 
absorption bands [Stiller et al, 1995], 

2.2. Data selection 

Data were screened prior to use to determine their quality, based on the following criteria: 
calibration against known standards such as the NOAA Climate Monitoring and Diagnostics 
Laboratory (CMDL) reference gases, the existence of intercomparison with other instruments 



Greenblatt et al. 


Tracer-based determination of vortex descent 


9 


[e.g., Chang et al., 1996; Herman et al, 1998; Michelsen et al, 1999; Toon et al, 1999; Moore et 
al, 2001], consistency of the N 2 0:CH 4 correlation with other instruments when both 
measurements were available for a single flight, and consistency of the profiles themselves with 
those of other instruments. Because data from dissimilar instruments were being compared, we 
have used total errors, consisting of precision and accuracy estimates, as well as uncertainties of 
calibration gases and/or spectroscopic parameters, in our analyses. 

Each data point of each flight was classified as vortex, extravortex, or mixed (and therefore 
excluded), using one of several categorization methods. For the balloon flights, each of which did 
not extend over large portions of the Arctic, hemispheric maps of modified Ertel’s potential 
vorticity (PV) [Lait, 1994], provided by the NASA Goddard DAO, were examined on several 9 
surfaces to determine vortex classification. To define the vortex edge, the Nash criterion [Nash et 
al, 1996] was used for flights from 1999-2000 and the Argus OMS in situ flight of June 30, 

1997; for the BONBON flight of November 30, 1991, the vortex edge was defined by the 
maximum PV gradient [Bauer et al, 1994]. 

For the unified N 2 0 data, flights were either classified as entirely extravortex from PV 
maps, or they were filtered to retain only vortex segments using the technique recently described 
by Greenblatt et al [2001], This technique used N 2 0 vs. 0 profiles measured inside the vortex to 
accurately determine the inner edge of the vortex boundary region; the inner edge was used to 
remove both large-scale flight segments of extravortex air, as well as small filaments of mixed 
composition. 

For the ASUR data, PV was converted to equivalent latitude (EqL) [Butchart and 
Remsberg, 1986] and points between 70°N and 90°N EqL were assumed to be inside the vortex. 
A comparison of the different months during 1999-2000 when the ASUR instrument was 
deployed showed that the vortex edge was located between 50 N and 70 N EqL. Each flight 
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used consisted of one or more profiles measured inside the vortex, from which a single weighted- 
average profile was constructed. 

HALOE data were filtered using a rotation/strain diagnostic called Q, which was 
integrated along rotational streamlines, with the vortex edge being defined as the rotational 
streamline with integrated Q of zero [Fair lie et al., 1999; Pierce et al., 2001], The data were 
further constrained to the high midlatitudes (> 45°N). Due to the small amount of data available 
at these latitudes, monthly average extravortex profiles were used for this paper. PV and Q were 
obtained from the United Kingdom Meterological Office (UKMO) UARS assimilated wind and 
temperature fields [Swinbank and O'Neill, 1 994], ' 

The ATMOS tracer data were sorted according to isentropic maps of PV and VMRs.of 
CH 4 , HNO 3 , 03 , and N 2 0. Averages included profiles with similar tracer distributions 
throughout the altitude range considered within selected geographic regions. Profiles displaying 
inversions (e.g., regions where 0 was not single-valued with respect to tracer VMR) were 
excluded from the averages. Classifications (vortex and midlatitude extravortex) were confirmed 
by comparison with tracentracer correlations [Michelsen et al, 1998a, 1998b]. 

2.3. Binning and inversion 

Filtered data were binned into 5 K 0 intervals in order to facilitate comparison among data 
sets. The tracer VMR reported for each 0 bin was calculated from the weighted mean of all data 
within the bin. Gaps in the data were filled in using linear interpolation. 

For the purpose of comparing changes in 0 on N 2 0 or CH 4 surfaces of constant VMR 
(isopleths), the data were “inverted” so that 0 could be expressed as a function of tracer VMR.. 
Inverted data were binned in intervals of 10 parts per billion by volume (ppbv) N?0 and 40 
ppbv CH 4 . The resulting tracer-binned data are more accurately described as “whisker” data, 
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since all data points whose VMR value ± 1 standard deviation (a) fell within a given tracer bin 
were considered in the calculation for that bin. The minimum and maximum 0 of all qualifying 
points were determined; 0 for the tracer-binned data point was defined as the midpoint of this 
range, and the la uncertainty in 9 was defined as half the range. Tracer-binned data were 
constructed from 0-binned data, rather than raw data, primarily to overcome data sparseness for 
some balloon flights. Gaps in the resulting data were filled in using linear interpolation. 

Figure 1 illustrates this idea, using a portion of the LACE N 2 O profile of November 19, 
1999 as an example. The region from 150-210 ppbv illustrates the need to consider the N 2 O 
uncertainties in the calculation of 0 for a given tracer bin, because N 2 O error bars in this region 
were much larger than the 10-ppbv bin size, and therefore single data points contributed to the 
calculation of 0 in several bins. Within the 230 and 240 ppbv N 2 O bins, where the 9-binned data 
display an inversion with respect to N 2 0, and at 300 ppbv and above, where the 0-binned data 
became essentially vertical, the midpoint of the range of 0 values within a bin served as the most 
useful definition of 0. The ±1 a envelopes of the N 2 0-binned data (dashed line) generally 
encompassed those of the 0-binned data (dotted line), which indicates that this approach 
provides a conservative estimate of the uncertainty relative to 0. The overestimate of the N 2 O- 
binned envelope with respect to the 0-binned envelope between 250-260 ppbv N 2 O was due to 
the lack of any 0-binned data points in these two N 2 O bins; the interpolation scheme smoothly 
filled in this region from surrounding data. 

2.4. Models 

SLIMCAT is a three-dimensional chemical transport model using horizontal winds and 
temperatures from the UKMO-UARS assimilation [Swinbank and O'Neill, 1994] to drive 
transport on an isentropic (0-based) vertical coordinate. Vertical motion is diagnosed by 
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calculating heating rates from a radiation model. The abundances of chemical species in the 
stratosphere are calculated using a detailed chemistry scheme [( Chipperfield , 1999]. The model 
simulates the distribution of all species involved in stratospheric ozone depletion, as well as 
some trace gases, such as N 2 0. The SLIMCAT model run started on October 21, 1991 and was 
integrated at low resolution (7.5' x 7.5' x 18 0 levels) until November 1999. Output from this 
run was then used to initialize a higher resolution simulation from November 3, 1999 to March 
31, 2000, on 24 0 levels ranging from 335-2726 K. It was assumed that an EqL greater than 70 'N 
was inside the vortex for the entire winter, and data between 70-80° N EqL were used to produce 
vortex averages. The N 2 0 data were recorded daily, but for the current analysis the data were 
averaged over 10 day periods, and interpolated as described above in 5 K 0 intervals. Data were 
also inverted into 10 ppbv N 2 0 bins. No error bars were calculated for the original data, so all 
data points were treated equally, and uncertainties in N 2 0 or 0 represented the observed 
variation for a given 0 or N 2 0 bin, respectively. 

REPROBUS is a three-dimensional chemical transport model whose advection is driven 
by the 6-hourly European Center for Medium-Range Weather Forecasting (ECMWF) 
meteorological analysis [ Lefevre et al., 1994, 1998]. The model includes a comprehensive 
treatment of gas-phase and heterogeneous stratospheric chemistry, as well as a semi-Lagrangian 
transport scheme integrated with 2° x 2° horizontal resolution, using all three components of the 
wind field. Calculations were performed on 42 vertical levels from the ground up to 0. 1 hPa. N 2 0 
and CH 4 tracer fields were initialized from zonal mean climatologies on November 1, 1999, and 
the model was run until March 30, 2000, recording N 2 0 and CH 4 once every 10 days. As with 
the SLIMCAT model, data between 70-80'N EqL were used to produce vortex averages, except 
for the two dates following the vortex breakup (March 20 and 30) when data were averaged 
between 80-90°N EqL, in order to obtain enough points within the vortex for good averaging 
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statistics. Average tracer VMR values were calculated on 5 K 0 intervals over the range 300-650 
K. Data were inverted into 10 ppbv N 2 O and 40 ppbv CH 4 intervals, the same as for the 
observational data. No error bars were originally calculated, so all data were treated equally, and 
uncertainties in 9 represented the observed variation for a given N 2 O or CH 4 bin. 

3. Results and Discussion 

3.1. Data from SOLVE/THESEO2000 

Figure 2 shows the balloon-borne data from 1999-2000. The data can be roughly divided 
into two groups: late fall, consisting of the November 19 LACE and December 3 MkIV and SLS 
profiles, and mid- to late winter, consisting of the January 27 and March 1 BONBON profiles 
and the March 5 LACE profile. The late fall profiles lie at higher 6 levels than the winter profiles 
over almost the entire range of N 2 O and CH 4 values, and is a result of descent. The late fall 
profiles also show a larger degree of variability than the later profiles, a point that will be 
addressed below. The March BONBON and LACE profiles, measured four days apart and 
launched from the same location (Esrange, Sweden), are in very good agreement with each other 
for both N 2 0 and CH 4 , with the exception of a difference seen in N 2 0, which falls significantly 
outside the la error bars beginning at 478 K. The cause of this discrepancy is not understood, 
but may be linked to a difference in the N 2 0:CH 4 correlations of these profiles below 40 ppbv 
N 2 0 (see below). The January 27 BONBON profile lies 9±6 K (la) above the two March 
profiles over 40-280 ppbv N 2 0 and 680-1600 ppbv CH 4 ; the difference is attributed to late 
winter descent. 

Figure 3 shows the unified N 2 0 data measured on the ER-2 in two panels, with a balloon 
profile shown for reference on each panel. Data are shown after binning into 5 K 0 intervals for 
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clarity. The profiles are very tight over the entire time period, with the ±lcr of N 2 O within a 
given 0 bin averaging 1 1 .4 ppbv over the range 370-470 K, the approximate vortex altitude region 
sampled by the ER-2. A decrease in 0 on a given N 2 O level can clearly be seen over the 52-day 
period, and often from flight to flight. 

Figure 4 shows ASUR profiles in three panels. A reference in situ balloon profile is also 
shown for each panel to indicate how closely the ASUR data agree with in situ measurements 
made at comparable times. A systematic overestimation of N 2 O at low VMR is apparent, which 
is due to the inherent limited vertical resolution of ASUR (5-10 km) and the high curvature of the 
profile in that region. At higher N 2 O values (e.g., above ~180 ppbv), where the change in slope is 
more gradual, the agreement with in situ measurements is much better. We have tried applying the 
ASUR averaging kernel to the balloon and unified N 2 O data, producing profiles with the same 
vertical resolution as ASUR. The resulting smoothed data showed much better agreement with 
ASUR at all N 2 O levels. Although using such an approach would allow comparison below 150 
ppbv, it was desirable to retain the higher vertical resolution of the other instruments. Therefore, 
in the analyses which follow, only the ASUR data on the 200 and 250 ppbv N 2 0 isopleths have 
been used. 

3.2. N20:CH4 correlation 

The N20:CH4 correlations for balloon data from 1999-2000 are presented in Figure 5. 

The five datasets show very good agreement with each other, except in the region between ~10- 
40 ppbv N 2 O, where the LACE flight of March 5 differs significantly from the other data, 
including the two BONBON flights. The difference in the LACE correlation may indicate real 
differences in the N20:CH4 correlation curve, because the MkIV profile of March 15 (not 
shown) closely follows the LACE correlation in this region. However, N 2 O values below 50 
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ppbv were not utilized in the work presented in this paper, so the disagreement does not impact 
our conclusions. 

Also shown in the Figure are the ATMOS vortex and midlatitude correlations [ Michelsen 
et al., 1998a], constructed from measurements made from 1992-1994. The data from 1999-2000 
fall in between the two ATMOS correlations between — 30- 1 50 ppbv N 2 O, and fall at higher CH 4 
values above —170 ppbv N 2 O. The secular trends in N 2 O and CH 4 [JPCC, 1996] are responsible 
for most of the differences at high N 2 O (i.e., lower altitudes), but at lower N 2 O values the 
influence of tropospheric trends is expected to be smaller because of the shorter photochemical 
lifetimes N 2 0 and CH 4 . The lack of sensitivity to such trends at these altitudes is supported by 
the excellent agreement between ATMOS correlations from 1994 and MkIV and ALIAS II 
correlations from 1997 [Hermann et al., 1998]. The differences shown in Figure 5 are consistent 
with differences between the ATMOS early fall developing vortex (or protovortex) and 
extra vortex correlations reported by Manney et al. [2000]. Such a difference can be attributed to 
some descent followed by mixing of extravortex air into the vortex, which is expected during this 
time of year [Manney et al., 2000; Plumb et al., 2000], The differences, however, are not as 
substantial as those between the extravortex correlation and the correlation observed in the spring 
vortex in 1993 (the ATMOS vortex curve in Figure 5). These differences have similarly been 
attributed to significant unmixed descent followed by mixing with extravortex air [Michelsen et al. 
[1998a, 1998b, 1999], 

For the purposes of this paper, it was not necessary to understand why the N 20 :CH 4 
correlation for the 1999-2000 winter differed from the ATMOS correlations, because the 
correlation was only needed to compare descent on N 2 O and CH 4 isopleths. Thus, a piecewise, 
weighted least-squares fit was generated, based on the five balloon flights, between 700-1800 
ppbv CH 4 (44-323 ppbv N 2 O). This fit is shown in Figure 5, with its functional form listed in 
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the caption. The fit was used to generate the CH 4 isopleths corresponding to the N 2 O isopleths 
displayed in Figure 6 . 

3.3. Late fall vortex inhomogeneity 

It is clear from examination of the three fall balloon profiles (Figure 2) and the fall ASUR 
data (Figure 4, panel (a)) that the vortex tracer profiles were quite variable in late fall. For 
instance, the mean MkIV-LACE difference was 23±3 K between 110-200 ppbv N 2 0, and 20±4 
K between 960-1200 ppbv CH 4 . The mean SLS-LACE difference was about half as large, 1 1±1 
K between 110-150 ppbv N 2 0. As MkIV and SLS were flown on the same balloon payload, the 
differences in their measurements were likely due to the fact that the instruments observed 
airmasses in different directions, corresponding to a horizontal separation of -800 km. 

The ASUR instrument sampled large regions of the vortex between November 30 and 
December 16. For this portion of the year, the data reveal a variability, indicated by 0±la, which 
was similar to the MkIV-LACE differences: 19 K at 200 ppbv N 2 0 and 25 K at 250 ppbv. For 
late winter ASUR profiles (February 27-March 15), the variability dropped to 17 K at 200 ppbv 
and 1 1 K at 250 ppbv (excluding the flight of March 8 , which exhibited unusually high N?0 
deviations; see magenta curve in Figure 4, panel (c)). By comparison, the unified N 2 0 variability 
in late winter (February 26-March 12) was 7.4 K at 200 ppbv and 9.8 K at 250 ppbv. Thus, 
though the noise in the ASUR data is certainly greater than for unified N 2 0, particularly at 200 
ppbv N 2 0, the large decrease in the ASUR variability at 250 ppbv is indicative of a more 
homogenous vortex in late winter than late fall. 

Further evidence for an inhomogenous vortex is given in a study by Ray et al. [2001], 

They found that they could explain the differences in tracer-tracer relationships between the 
November and March LACE profiles with a simple differential diabatic descent model. In this 
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model, air parcels underwent different amounts of descent prior to the November 1 9 
measurement, resulting in air with significantly different tracer VMRs on the same 0 surfaces, 
which subsequently mixed [see e.g., Plumb et al., 2000]. Their model simultaneously fit six 
measured stratospheric tracers, including N 2 O and CH 4 , producing a self-consistent result with 
differential diabatic descent ranging from ~2 km at 1 5 km to ~6 km above 30 km. Thus, their 
analysis indicated that the vortex was inhomogenous in the fall but became well mixed later in the 
winter. 

Manney et al. [ 2000 ] cited evidence for vortex inhomogeneity even earlier in the winter of 
1994-1995. From November 3-12, 1994, the ATMOS instrument observed several laminae 
within the developing polar vortex whose ^OiCFLj correlations indicated both entrainment of 
tropical air, as well as midlatitude air that had undergone significant descent. The tracer vs. 0 
profiles indicated large deviations from the bulk correlation in these regions, resulting in a vortex 
with substantial inhomogeneity. 

All of the above indicates that the vortex tends to be fairly inhomogenous in late fall. By 
early November, the vortex is typically well established only above 600 K [Waugh and Randel, 

1999] , and in 1999 had only recently formed in the lower stratosphere [Manney and Sabutis, 

2000] . Thus, it is a reasonable conclusion that complete mixing had not yet occurred. The smaller 
variability in 0 for a given tracer isopleth seen later in the winter supports the idea of a more 
uniformly mixed vortex. 

3.4. Late fall average profile 

For the purpose of estimating an average late fall vortex profile with realistic 
uncertainties, we have combined the LACE and MkIV balloon profiles in the following manner. 
(The SLS profile was not included because of its large error bars at low and high N 2 O levels. 



Greenblatt et al. 


Tracer-based determination of vortex descent 


18 


However, the SLS measurements were bounded by the other two profiles, so their omission did 
not detract from the goal of creating a profile that spanned the range of observations. The ASUR 
data were not included due to their low vertical resolution). It was not appropriate to calculate 
error-weighted averages because the profiles did not look at the same air parcels. Thus, for each 0 
bin, the minimum and maximum tracer VMR ± la errors were determined; the “average” VMR 
was defined as the midpoint of that range, and the la uncertainty was defined as half the range. 
The procedure was applied independently to the N 2 0 and CH 4 data. The date assigned to each 
profile was the average date of the two datasets, November 26, 1999. The profiles are shown in 
Figure 7. 

3.5. Time series 

The balloon, unified N 2 0 and ASUR data have been combined and presented as a time 
series along several N 2 0 and CH 4 isopleths in Figure 6. Data shown here have been inverted as 
described in the previous section to generate 9 as a function of N 2 0 or CH 4 . The levels of 50, 

100, 150, 200, and 250 ppbv N 2 0 were chosen to span the useful range of the measurements: at 
levels below 50 ppbv, there is almost no unified N 2 0 data, and the balloon profiles become quite 
steep in 0, making comparisons along isopleths difficult; much above 250 ppbv, the curves 
cluster on top of one another at all time periods, indicating marginal descent and probably 
significant mixing with midlatitude air on these isopleths [e.g., Abrams et al., 1996a], CH 4 levels 
of 720, 920, 1080, 1280, and 1480 ppbv were chosen to correspond approximately (to within 
one 40-ppbv CH 4 bin) to these N 2 0 levels, based on the average N 2 0:CH 4 correlation derived 
from the balloon measurements, shown in Figure 5. Simple fits through the data (linear least- 
squares fits in dense data regions (November-December, January-February, and February-March 
deployments), and connected points in between these regions) have been applied to the N 2 0 data 
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to indicate the general progress of the descent throughout the winter. These same fits have been 
reproduced on the CH 4 panel for reference. The more limited CH 4 data are very consistent with 
these fits throughout the winter. 

The progress of descent is evident in detail in Figure 6. The rates of descent generally 
increased with decreasing tracer VMR, as indicated by the slope of the fits, and the greatest rates 
of descent occured before January 1. Between the midwinter (January-February) and late winter 
(February-March) data segments, the descent rates slowed considerably. Data from all 
instruments were in excellent agreement with each other throughout the winter, with the 
exception of the LACE and MkIV late fall data, which lay outside each other’s error bars on the 
150 ppbv N 2 0 and 1080 ppbv CH 4 isopleths because of the vortex inhomogeneity discussed 

earlier. 

3.6. Vortex starting profile 

Many studies of vortex descent have calculated the total descent over the winter season 
starting from profiles measured in late fall [e.g., Bauer et al., 1994; Strahan et al., 1994; Muller et 
al., 1996; Hartmann et al., 1997]. However, it is known that descent can begin quite early in the 
winter season. For instance, Michelsen et al. [1998b] and Manney et al. [1999] showed that 
profiles measured inside the Arctic vortex during November 3-12, 1994 had already descended 
several km relative to extravortex profiles. Abrams et al. [1996a] calculated from analysis of PV 
contours that descent had begun in the fall of 1992 in the Arctic between October 1 and 
November 12, with the earliest dates occurring at the highest altitudes (> 42 km). Rosenfield et al. 
[1994] chose a starting date of November 1 for model simulations of Arctic vortex descent in 
1988-1989 and 1991-1992, with significant descent rates calculated from the first day of the 


simulations. 
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By November 1 9, 1 999, the date of the first LACE profile, significant descent had 
occurred at all altitudes relative to extravortex profiles, as will be demonstrated shortly. Thus, 
using this profile (or the November 26, 1999 average fall profile) as the starting point for 
calculating descent would miss a significant amount of descent that occurred earlier in the fall. In 
order to determine the total vortex descent in 1999-2000, a reasonable starting profile had to be 
established. Prior to vortex formation in late fall, air in polar regions is constantly mixed over 
horizontal distances of several thousand km [Schoeberl and Hartmann, 1991], Thus, the tracer 
vr. 0 profile of polar air before vortex formation strongly resembles air of midlatitude 
composition. This has been confirmed observationally by a number of studies [Bauer et al., 1994; 
Abrams et al., 1996a; Herman et al., 1998; Jost et al., 1998; Michelsen et al., 1998a], It has been 
recognized that the tracer vs. 0 profiles exhibit small seasonal changes [ Abrams et al., 1 996a], and 
that tropical profiles appear significantly different due to their partial isolation from the 
midlatitudes [e.g., Jost et al., 1998; Herman et al., 1998; Michelsen et al., 1998a], However, by 
limiting our examination to high-latitude, extravortex regions, we have assembled a selection of 
candidate vortex starting profiles. 

These profiles are shown in Figure 7, along with the November 26, 1999 average vortex 
profile for reference (solid black line). The extravortex profile from the BONBON instrument, 
shown as a solid red line, was launched from Kiruna, Sweden on November 30, 1991 [Bauer et 
al, 1994], Averaged extravortex profiles from the April 1993 and November 1994 ATMOS 
missions are shown as orange and green solid lines, respectively (H. A. Michelsen, private 
communication). The Argus profile, shown as a solid cyan line (for N 2 0 only), was measured in 
Fairbanks, Alaska on June 30, 1997 [Jost et al., 1998]. It exhibited several large excursions to low 
N 2 0 that have been interpreted as vortex remnant filaments, and were therefore omitted from 
Figure 7. The limited extravortex data from 1999-2000 included a unified N 2 0 profile from 
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March 16 (dark blue solid line, N 2 O only), and monthly averaged HALOE profiles from 
December, January, and March poleward of 45’N (magenta, black, and red, respectively, with 
dotted lines; CH 4 only). 

While there is variation among profiles of up to ~70 ppbv N 2 O and —400 ppbv CH 4 on 
certain 0 levels, the overall pattern is consistent through multiple years and seasons, and is, most 
importantly, quite distinct from the November 26, 1999 average vortex profile. Therefore, we 
assert that any profile from the set presented here can be used as a proxy for the vortex starting 
condition for 1999-2000. The November 1994 ATMOS average was chosen, which when both 
N 2 O and CH 4 were considered, occupied a roughly median value at all 0 values. We denoted this 
profile as the “START” profile, and used it in all subsequent calculations. 

3.7. Descent calculations 

In order to determine the amount of descent between measurements made at different 
times, differences in 0 between profiles have been calculated, with the results shown in Figure 8 . 
Vertical bars indicate ±1 ct uncertainties, which were calculated as the quadrature sum of 
individual profile uncertainties. The diamonds/solid lines represent differences (A0, expressed in 
K) between the START profile and the November 26, 1999 average; squares/dotted lines, 
between the November 26 average and the January 27, 2000 BONBON profile; and 
triangles/dashed lines, between the BONBON profile and the March 5, 2000 LACE profile. 
Differences are shown for both N 2 O and CH 4 data in each case. 

For the START-November 26 difference, A0 ranged from 397±15 K (la) at the lowest 
VMR values (30 ppbv N 2 O and 640 ppbv CH 4 ), to 28±13 K above 200 ppbv N 2 O and 1280 
ppbv CH 4 . These differences were larger than those observed in any other profile pair, and 
indicated that a great amount of descent had occurred by late fall 1 999. Such large descent was 
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also inferred for early November 1994, when the ATMOS instrument examined air both inside 
and outside the forming polar vortex [Michelsen et al., 1998b]. While that study expressed 
descent in terms of changes in r only, differences in terms of 0 have been calculated and found to 
be comparable to our results: 220-240 K at 30 ppbv N 2 O and 640 ppbv CH 4 , and 380 K or more 
below 10 ppbv N 2 O and 390 ppbv CH 4 (H. A. Michelsen, private communication). 

The November 26-January 27 differences were much smaller than the START-November 
26 differences below 170 ppbv N 2 O and 1160 ppbv CH 4 , but were fairly comparable above 
these isopleths. From 40-280 ppbv N 2 O and 680-1600 ppbv CH 4 , which covered most of the 
tracers’ dynamic ranges, A0 was fairly constant at 54±10 K. Between January 27 and March 5, 
A0 was smaller still, only 10±6 K over the same ranges. The differences between the START 
profile and the March 5 LACE profile (not shown Figure 8 ), as an indication of the total descent 
during 1999-2000, ranged from 501±21 K at 30 ppbv N 2 O and 640 ppbv CH 4 to 49±23 K above 
200 ppbv N 2 O and 1280 ppbv CH 4 . 

It is important to emphasize that these descent results represent first-order estimates 
only, because the effects of mixing have not been accounted for. “Mixing” here indicates the 
process of vortex homogenization, incorporating air that has either become entrained from 
midlatitudes, or undergone differential diabatic descent within the vortex. In either situation, 
adjacent air masses will have undergone different amounts of descent, and tend to mix along 0 
surfaces. Both of these processes can be active during vortex formation in the fall, and can 
substantially alter tracentracer as well as tracer:© relationships [Plumb et al., 2000], The general 
effect is to reduce the apparent descent [Abrams et al., 1996b; Muller et al, 1996]. However, 
once the vortex has become isolated from the rest of the stratosphere, which is typically 
complete by late fall, the “no mixing” assumption is a fairly good one. For instance, in the study 
by Ray et al. [2001], which explicitly modeled mixing as a result of differential diabatic descent in 
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the early fall, only small changes (< 20 K) in the calculated descent between the November 19, 
1999 and March 5, 2000 LACE measurements were seen relative to the simple case where mixing 
was ignored. 

Descent rates are shown in Figure 9, panels (a)-(e), as well as in Table 2. Rates are shown 
for all profile pairs shown in Figure 8, with the exception of the START-November 26 
difference, since the effective 1999 date of the START profile was not known. Rates for several 
unified N 2 O profile pairs are also shown. Although the unified N20-based rates display generally 
higher error bars than those based on the balloon data, the latter were included because their 
higher time density enabled a more detailed temporal reconstruction of the descent rates. There is 
a general decrease in the descent rate with time for all N 2 O isopleths shown, ranging from 
between 0.56±0.18 to 1.05±0.21 K/day in the earliest time interval (November 26-January 27), 
to -0.17±0.37 to 0.33±0.47 K day in the latest interval (February 26-March 12). This is entirely 
expected, and is a reflection of increasing solar radiation slowly reducing the net radiative cooling 
as winter passes into spring. There is also a general decrease in the rate with increasing N 2 O 
isopleth: e.g., for the time interval January 27-March 5, the rate decreases from 0.39±0.15 K/day 
at 50 ppbv to 0.13±0.19 K/day at 250 ppbv. This N 2 O trend reflects the greater radiative cooling 
at higher altitudes, while the much lower descent by 250 ppbv is indicative of the increased 
mixing near the tropopause. However, a notable exception is the 200 ppbv isopleth, which has 
the highest descent rate in the earliest interval. 

Weighted-average descent rates across N 2 O isopleths for each time interval are also 
shown in Figure 9 panel (f), as well as in Table 2, decreasing from 0.82±0.20 K/day between 
November 26-January 27 to 0.10±0.25 from February 26-March 12. Finally, a winter-average 
descent rate (from November 26-March 5; results for November 26-March 12 are very similar) is 
shown in Table 2 for each N 2 O isopleth, as well as a grand average. A decrease from 0.75±0.10 
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K/day at 50 ppbv to 0.40±0.1 1 K/day at 250 ppbv is observed, with the average winter 1999- 
2000 descent rate across 50-250 ppbv N 2 O at 0.61±0.13 K/day. 

3.8. Comparison with other studies 

Descent rates have been reported previously relative to 9 [Schoeberl et al., 1 992; 
Rosenfield et al, 1994; Strahan et al., 1994; Lucic etal., 1999], altitude (z) [ Schoeberl et al., 1992; 
Rosenfield et al., 1994; Traub et al., 1995; Abrams et al., 1996a], and/or tracer isopleth [Bauer et 
al, 1994; Muller et al., 1996; Hartmann et al., 1997], and over a wide range of time periods. In 
order to compare our results with previously reported values, we have related values found in the 
literature to N 2 O isopleths, as summarized in Table 3. The methods used for each study are 
explained in the footnotes to Table 3. 

Descent rates originally reported relative to 0 are generally consistent with our values (see 
Table 2). Strahan et al. [1994] and Muller et al. [1996] reported observation-based net winter 
(e.g., December-March) descent rates for 1991-1992 and 1991-1992 through 1995-1996, 
respectively, between 0.5-0. 7 K/day over 100-230 ppbv N 2 O. The 1991-1992 model study by 
Lucic et al. [1999] reported net winter descent rates of 0.55-0.9 K/day between 70-250 ppbv 
N 2 O. All three studies were in line with the 1999-2000 winter averages from November 26- 
March 5 (see Table 2). Schoeberl et al. [1992] reported both observed and modeled descent rates 
for a ~1 -month period starting in early January 1992. Their results, 0.09±0.16 to 0.36±0.05 
K/day between 1 10-240 ppbv N 2 O, were well within the scatter of calculated rates for the 
comparable period January 23-February 2, 2000. Strahan et al. [1994] reported modeled monthly 
descent rates from October through March (see Table 3); the rates at 150 ppbv N 2 O from 
December through March (0.2-0. 9 K/day) are in good agreement with the 1999-2000 results 
(0.33±0.47 to 0.69±0.25 K/day). Rosenfield et al. [1994] similarly reported modeled descent rates 
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over most of the 1988-1989 and 1991-1992 winter seasons (November 1-March 21); they are 
reported for two time periods in Table 3. While the rates for November 1 -January 22 across 50- 
250 ppbv N2O (0.72-1.20 K/day) compare well with the 1999-2000 rates for November 26- 
January 27 (0.56±0.18 to 1.05±0.21 K/day), the rates for January 22-March 21 (0.60-0.95 
K/day) are significantly higher than the similar 1999-2000 time period of January 27-March 5 
(0.07±0.09 to 0.39±0.15 K/day). 

Bauer et al. [1994] reported observation-based descent rates from October 31, 1991 to 
January 18, 1992 of 2. 7-3.4 K/day at 50-100 ppbv N2O, which are far higher than those 
calculated for the similar period November 26, 1999-January 27, 2000 (0.81±0.21 to 0.97±0.18 
K/day). The difference may be partly due to the earlier starting time of the calculation, 
interannual variation in descent rates, or both. Strahan et al. [1994] reported a model descent rate 
for October 1991 of 1.2 K/day at 250 ppbv N2O; although not reported in the study, the rate at 
50-100 ppbv would probably be much larger, consistent with the trend in that paper of higher 
descent rates on smaller N2O isopleths. However, Rosenfield et al. [1994], who modeled the 
descent for 1988-1989 and 1991-1992, began their simulation on November 1, almost coincident 
with the first Bauer et al. [1994] profile, and reported substantially smaller rates though January 
22 (1 .02-1.20 K/day at 50-100 ppbv N2O). While no systematic, multi-year analysis of Arctic 
descent rate variability has been performed, a study of Antarctic HALOE data from 1992-1997 
by Kawamoto and Shiotani [2000] revealed interannual changes in descent rates of ±20% (39-59 
m/day at 600 ppbv CH4). This amount of variability is not nearly enough to account for the 
differences from the 1999-2000 results, however. 

It was hoped that by examining studies that reported descent rates in terms of both 0 and 
z, a simple correspondence could be derived between the two sets of rates. However, the data 
showed little or no relationship between the two coordinate systems; for instance, Rosenfield et 


Greenblatt et al. 


Tracer-based determination of vortex descent 


26 


al. [1994] reported a substantial increase in the z-based descent rates after January 22, whereas 
the 0-based rates in the same study decreased markedly. Another example is Bauer et al [1994], 
whose z-based rates, while on the high end of the range reported by other studies, did not stand 
out nearly as much as the 0-based rates discussed earlier (see Table 3). Thus, the relationship 
between 0 and z appears complex, and we did not attempt to make any comparison between 
studies that reported descent rates in terms of z and the 0-based 1999-2000 descent rates. 

3.9. Model comparisons 

Figure 10 shows the comparison between the REPROBUS model and fits of the 
observational data. The agreement differed considerably from N 2 O to CH 4 , with each tracer 
exhibiting better agreement on some isopleths and worse agreement on others. In terms of 
absolute differences, the N 2 O model data agreed most favorably with the observational fits from 
150-250 ppbv, while the CH 4 data agreed quite well for the 720 ppbv isopleth (equivalent to 50 
ppbv N 2 O), but had worse agreement for all other isopleths. Most of the differences observed 
may be due to poor model tracer initializations and/or inadequate model tracer photochemistry. 
However, while the differences between model and observation may be considerable for several 
isopleths, they did not change significantly with time (with the exception of the 50 ppbv N 2 O 
and 1480 ppbv CH 4 isopleths, whose model-observation differences increased markedly with 
time). This relative agreement indicates that the model-derived descent rates were, in fact, fairly 
accurate throughout the winter. The disagreement at 50 ppbv N 2 O may indicate that descent was 
not strong enough in the model at these altitude levels, while the disagreement at 1480 ppbv CH 4 
was probably due to the close proximity to the bottom of the vortex, which is typically located 
near 400 K [e.g., Manney et al., 1 994], 

A comparison of the SLIMCAT N 2 O results and fits of the observational data is shown 
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in Figure 11. Offsets between model and observation were of the same magnitude (~20 K) as the 
REPROBUS differences for the 100-250 ppbv isopleths, and the 50 ppbv isopleth agreed much 
better with the observed data than did the 50 ppbv N 2 0 REPROBUS isopleth. There was a 
slight change in the offsets between 100-200 ppbv around day 27, resulting in better agreement 
between model and observation, while the 250 ppbv isopleth offset increased slightly over time. 
The 50 ppbv isopleth, despite its larger error bars, appeared to track the observed data best, with 
the exception of days 37 and 87, where error bars were much higher; these increases were 
associated with a more distorted vortex, resulting in a significantly wider range of VMRs sampled 
in the 70-80* N EqL region. Despite these 0 offsets, the good agreement in slopes between the 
SLIMCAT results and the observations indicated that the descent rates calculated were very 
similar to those derived from our study. 

4. Summary 

The descent of the 1999-2000 Arctic vortex has been quantified using N 2 0 and CH 4 
tracer measurements from a number of instruments. Coverage spanned the last stages of vortex 
formation in late November 1999 through vortex breakup in mid-March 2000, with two large 
gaps (December 17-January 13, and February 4-25). Changes in 0 were determined on five N 2 0 
and CH 4 isopleths, and descent rates were calculated on each N?0 isopleth for several time 
intervals throughout the fall and winter. From November 26-March 5 the vortex descended 
65±12 K when averaged over 40-280 ppbv N 2 0 and 680-1600 ppbv CH 4 . The maximum descent 
rates occurred in the late fall/early winter phase: between November 26 and January 27, the 
calculated rate was 0.82±0.20 K/day averaged over 50-250 ppbv N 2 0. By late winter (February 
26-March 12), the average rate had decreased to 0.10±0.25 K/day. Descent rates also decreased 
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with increasing N2O isopleth; for instance, the winter average (November 26-March 5) descent 
rate varied from 0.75±0. 10 K/day at 50 ppbv to 0.40±0.1 1 K/day at 250 ppbv. The winter 
average rate over 50-250 ppbv N2O was 0.61±0.13 K/day. 

Larger variability in 0 on tracer isopleths in late fall vs. winter indicated that the vortex 
was much less homogenous in the late fall. The degree of inhomogeneity, derived from several 
instruments (LACE, MkIV, SLS, and ASUR), was 19-23 K between 1 10-200 ppbv N 2 0 and 
960-1200 ppbv CH4. This inhomogeneity decreased significantly in the March datasets; the 
unified N2O variability averaged only 7.4 K at 200 ppbv N2O. Because of the vortex 
inhomogeneity in late fall, an average fall vortex profile was constructed from the November 19 
LACE and December 3 MkIV measurements, which encompassed the observed variability of all 
the fall measurements. This profile, referred to as the November 26, 1999 vortex average, was 
used in subsequent calculations of vortex descent. 

A comparison of high-latitude, extravortex profiles measured in different years - and 
seasons, while revealing significant VMR differences on a given 0 level, also showed that 
extravortex profiles were more similar to each other than to the early vortex profiles. Therefore, 
the median profile, an average extravortex profile measured by ATMOS in November 1994, was 
chosen as the vortex starting profile for 1999-2000 (denoted “START”) and was used to estimate 
the amount of descent prior to the first vortex measurements in November and December. The 
difference between the START profile and the November 26 average profile showed significant 
descent prior to the first profiles measured in fall 1999-2000, up to 397±15 K at the lowest N2O 
and CH 4 levels (30 ppbv N 2 0 and 640 ppbv CH 4 ), and falling to 28±13 K for N 2 0 larger than 
200 ppbv and CH 4 larger than 1280 ppbv. Such early descent has been inferred for prior years in 
the Arctic, e.g., 1988-1989 and 1991-1992 [Rosenfield et al., 1994], 1993-1994 [Abrams et al, 
1996a] and 1994-1995 [ Michelsen et al., 1998b]. While these differences did not include the 
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effects of mixing, as first-order approximations they indicated that a great deal of descent had 
occurred prior to the late fall vortex measurements. 

Descent rates from previous 9-based studies corroborate the results of our study very 

well across a wide range of time periods and N 2 O isopleths, with the exception of one 

observational study (Bauer et al. [1994]) and the late winter time period of one model study 

(. Rosenfield et al. [1994]). Because it has been shown, both with our own results and those of 

prior studies, that descent rates increase substantially with decreasing N 2 O or increasing altitude, 

it is more meaningful and useful to calculate rates for several tracer isopleths (as we have done), 9 

levels, or z levels, than to report a single number. Likewise, because of the large change in descent 
¥ 

rates observed throughout the fall and winter, a single, season-average descent rate is not as 
informative as rates reported for several selected times. 

Both the REPROBUS and SLIMCAT models calculated descent rates that agreed quite 
favorably with observations. Although both models showed disagreements with observations of 
up to ~20 K (except for the 50 ppbv N 2 O isopleth of REPROBUS), the causes were somewhat 
different for each model. For REPROBUS, the quality of the tracer initializations directly 
affected the agreement, as there was no “spin up” period in the model. For SLIMCAT, where 
tracer distributions were initiated in 1991, the differences with observation were due almost 
exclusively to model chemistry and transport, which over the eight year run time of the model 
should be regarded as very good. Given the different approaches employed by the two models 
for deriving these rates (direct use of the vertical wind field component for REPROBUS vs. a 
radiative model for SLIMCAT), the success of both models in reproducing the observed descent 
via trace gas measurements should be highlighted. 
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Figure Captions 

Figure 1. Illustration of the inversion technique, using part of the N 2 0 profile from LACE on 
November 19, 1999 as an example. Diamonds represent “0-binned” data, that is, data 
which have been binned into 5 K 0 intervals, with horizontal bars indicating ±1 ct 
uncertainties in N 2 0 for each data point. Squares represent inverted “N 2 0-binned” data, 
which have been placed into 10 ppbv N 2 0 intervals from the 0-binned data, with vertical 
bars indicating ±1 a uncertainties in 0 for each data point. The ±la envelopes of the N 2 0- 
binned data (dashed line) generally encompassed those of the 0-binned data (dotted line), 
which indicates that this approach provides a conservative estimate of the uncertainty 
relative to 0. The overestimate of the N 2 0-binned envelope with respect to the 0-binned 
envelope between 250-260 ppbv N 2 0 was due to a lack of any 0-binned data points in 
these two N 2 0 bins; the interpolation scheme smoothly filled in this region from 
surrounding data. 

Figure 2. Balloon profiles from 1999-2000 for (a) N 2 0 and (b) CH 4 . The legend indicates the date 
of each flight (YYYYMMDD format) and the instrument acronym. The color scheme is 
the same for both N 2 0 and CH 4 . Note that SLS only measured N 2 0. Horizontal bars 
indicate ±1 a uncertainties of each data point. Other lines connect the points of each 
profile. 

Figure 3. Unified N 2 0 (“uN 2 0”) profiles from 1999-2000 for the two ER-2 deployments: (a) 

January 20-February 3, 2000 and (b) February 26-March 12, 2000. In addition to unified 
N 2 0 data, the LACE balloon profile of March 5, 2000 (see Figure 2) is included in each 
panel to facilitate comparison. All data have been binned into 5 K 0 intervals for clarity. 
Horizontal bars indicate ±1(7 uncertainties, and other lines connect the points of each 
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profile. As the ER-2 aircraft, from which unified N 2 O data was measured, had a flight 
ceiling of ~21 km, the data do not extend as high in 0 as for the balloon measurements. 

Figure 4. ASUR profiles from 1999-2000 for the three DC -8 deployments: (a) November 30- 
December 16, 1999, (b) January 14-29, 2000 and (c) February 27-March 15, 2000. In 
addition to ASUR data, one balloon flight from Figure 2 is included in each panel to 
facilitate comparison. Florizontal bars indicate ±la uncertainties, and other lines connect 
the points of each profile. The disagreement between the ASUR and balloon data for N 2 O 
< -150 ppbv is due to the limited vertical resolution (-5-10 km) of the ASUR instrument, 
and the high curvature of the profiles in this region. 

Figure 5. N 2 0 :CFl 4 correlations for balloon flights from 1999-2000. Symbols indicate individual 
flights. Horizontal and vertical bars indicate ±la uncertainties. Shown also are ATMOS 
correlations for midlatitude and polar vortex air [Michelsen et al., 1998a], ^CEGFLt fit 
function is indicated by heavy solid line, with equation [N 2 O] =po + /?i [CH 4 ] -fc 
/? 2 [CH 4 ] 2 , where [N 2 O] and [CH 4 ] are in ppbv. Coefficients p-, are defined piecewise for 
two CH 4 regions: [CH 4 ]=700-1300 ppbv, p=[-\ 03.265, 0.172407, 5.32430x10-5]; 
[CH 4 ]= 1300- 1800 ppbv, />=[-80.07 1 5, 0.223859, 0], 

Figure 6 . Descent time sequence for five N 2 O (a) and CH 4 (b) isopleths. Each isopleth is denoted 
by a separate color. CH 4 isopleths correspond to N 2 O isopleths by the equation in Figure 
5. Symbols indicate individual instrument measurements. Vertical bars indicate ±la 
uncertainties. Dotted lines are fits for each N 2 O isopleth (fits are repeated on CH 4 panel), 
where linear least-squares fits were made in dense data regions (November-December, 
January-February and February-March deployments), and lines are connected together to 
span these regions. “Fall vortex” profiles (crosses) are constructed from averages of the 
November 19 LACE and December 3 MkIV profiles (see text), which are shown in their 
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entirety on Figure 7. 

Figure 7. High-latitude, extravortex profiles from several years and seasons for (a) N2O and (b) 
CH4. Also shown for reference is the November 26, 1999 “fall vortex” profile, whose 
construction is explained in the text. The ATMOS 1994 average profile occupies 
approximately a central location among extravortex profiles, when both N2O and CH4 are 
considered. This profile, denoted “START,” was used as a proxy for the vortex starting 
profile in subsequent calculations. 

Figure 8. Differences between profiles for (a) N 2 0 and (b) CH 4 , indicating descent (assuming no 
mixing). Vertical bars indicate ±lo uncertainties, which are calculated as the quadrature 
sum of individual profile uncertainties. Legend: START — November 26, 1999 fall vortex 
(diamonds/solid line), November 26, 1999 fall vortex - January 27, 2000 BONBON 
(squares/dotted line), January 27, 2000 BONBON — March 5, 2000 LACE 
(triangles/dashed line), zero line (dotted horizontal line). 

Figure 9. Descent rates (K/day) on different N2O isopleths for pairs of profiles on five N2O 

isopleths in panels (a)-(e), and average descent rates in panel (f). Horizontal bars indicate 
time differences, and vertical bars indicate ± 1 c uncertainties. Legend: November 26, 1999 
fall vortex - January 27, 2000 BONBON (diamond/solid line), January 27, 2000 
BONBON - March 5, 2000 LACE (triangle/solid line), January 20-February 2, 2000 
unified N2O (square/dashed line), February 2-February 26, 2000 unified N2O (X/dashed 
line), February 26-March 12, 2000 unified N 2 0 (diamond/dashed line), zero line (dotted 
horizontal line). 

Figure 10. Comparison of fits of observed data with the REPROBUS model for (a) N 2 0 and (b) 
CH 4 . Colors denote different isopleths following the same scheme as in Figure 6. Symbols 
indicate 1-day snapshots, with vertical bars indicating ±1 g uncertainties. Dotted lines are 
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the N 2 O fits from Figure 6. 

Figure 1 1 . Comparison of fits of observed data with the SLIMCAT model for N 2 O. Colors 

denote different isopleths following the same scheme as in Figure 6. Symbols indicate 10- 
day averages, with vertical bars indicating ±lcr uncertainties. Dotted lines are the N 2 O fits 
from Figure 6. 


Tables 

Table 1. Instruments used in this study. 

Data utilized 

Instrument Platform * Location ! N 2 0 CH £ 0 Patels! tYYYY.MM.DD> 

A. Measurements from SOLVE/THESE02000 


LACE 

OMS in situ 3 

Esrange 

Y 

Y 


1999.11.19, 2000.03.05 

JPL Ozone 

OMS in situ* 

Esrange 



Y 

1999.11.19 • 

NOAA Hygrometer OMS in situ 3 

Esrange 



Y 

2000.03.05 

BONBON 

TRIPLE 3 

Esrange 

Y 

Y 

Y 

2000.01.27, 2000.03.01 

Argus* 

ER-2 b 

Kiruna 

Y 



2000.01.20-2000.03.16 

ALIAS* 

ER-2b 

Kimna 

Y 



2000.01.20-2000.03.16 

ACATS* 

ER-2b 

Kiruna 

Y 



2000.01.20-2000.03.16 

WAS* 

ER-2b 

Kiruna 

Y 



2000.01.20-2000.03.16 

MMS 

ER-2b 

Kiruna 



Y 

2000.01.20-2000.03.16 

MkIV 

OMS remote 3 

Esrange 

Y 

Y 

Y 

1999.12.03 

SLS 

OMS remote 3 

Esrange 

Y 


Y 

1999.12.03 

ASUR 

DC-8b 

Kiruna 

Y 


Y 

1999.11.30-2000.03.15 

HALOE 

UARSc 



Y 

Y 

1999.10.01-2000.03.31 
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B. Additional measurements from prior years 
BONBON CNESa Esrange Y 

ATMOS Space Shuttle Y 

Argus OMS in situ a Fairbanks Y 

NOAA Hygrometer OMS in situ a Fairbanks 
*Platform types: a Balloon ^Aircraft ^Satellite 
fBase locations: Esrange, Sweden (68”N, 20°E); Kiruna, Sweden (68°N, 21°E); Fairbanks. Alaska 
(65 ’N, 148'W). 

JData from the four N 2 O instruments on board the ER-2 (Argus, ALIAS, ACATS and WAS) 
were combined into a unified N 2 O data product; see Hurst et al. [2001]. 

Table 2. Calculated descent rates (K/day) for the 1999-2000 winter. Uncertainties are la. 

N 2 O isopleth (ppbv) 


Ranee of dates 

50 

100 

150 

200 

250 

Average 

A. Balloon measurements 






November 26-January 27 

0.97±0.18 

0.81±0.21 

0.69±0.25 

1.05±0.21 

0.56±0.18 

O 

bo 

to 

H- 

p 

k> 

0 

January 27-March 5 

0.39±0. 15 

0.33±0.09 

0.33±0. 15 

0.07±0.09 

0.13±0.19 

0.24±0.14 

B. Unified N 2 O measurements 






January 23 -February 2 

n/a 

0.75±0.35 

0.00±0.35 

0.50±0.35 

0.25±0.79 

0.41±0.34 

February 2-February 26 

n/a 

0.31±0.15 

0.52±0.23 

0.31±0.23 

0.21±0.44 

0.35±0.16 

February 26-March 12 

n/a 

0.17±0.24 

0.33±0.47 

-0. 17±0.37 

0.00±0.53 

0.10±0.25 

C. Winter average 







November 26-March 5 

0.75±0.10 

0.63±0.13 

0.55±0.16 

0.68±0.13 

0.40±0. 1 1 

0.61±0.13 


Y Y 1991.11.30 

Y Y 1993.04.08-1993.04.16, 

1994.11.03-1994.11.12 
1997.06.30 
Y 1997.06.30 
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Table 3. Descent rate comparison with other studies. Uncertainties are la. Potential temperature 
is denoted by 0; altitude, by z. Numbers without footnotes indicate they were taken directly 
from the text or table in the reference. 







n 2 o 

Descent rate 

Reference 

Yeaifs) 

Time period 

OIK) 

z Ikin') 

(ppby) 

fK/davl 

Cm/dav) 

A. Observations 







Schoeberl 

1988-9 

Jan. 3-Feb. 10 

450a 

20a 

110b 

0.36±0.05 a 

98±20a 

et a/. [1992] 



420a 

18a 

160b 

0.27±0.11a 

69±33a 




390a 

16.5a 

210b 

0.09±0.16a 

12±35a 

Bauer 

1991-2 

Oct. 31 -Jan. 18 

510* 

21-0i 

50 

3.4= 

97= 

etal. [1994] 



445* 

18.4* 

100 

2.7= 

113=. 

Strahan 

1991-2 

Dec. -Feb. 

410 

n/a 

160-230 d 0.5 

n/a 

etal. [1994] 



470 

n/a 

1 00- 1 40 d 0.7 

n/a 

Traub 

1991-2 

Jan.-Feb. 

n/a 

18 

120 e 

n/a 

51±8 

etal. [1995] 








Abrams 

1992-3 

Oct./Nov. -Apr. 

n/a 

30 

55f 

n/a 

62 

etal. [1996a] 



n/a 

24 

120 f 

n/a 

39 




n/a 

20 

180 f 

n/a 

24 
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Muller 

1991-2 

Nov./Dec.-Apr. 

420> 

n/a 

150s 

0.5-0.6 h 

40h 

etal. [1996] 

to 1995-6 






Hartmann 

1994-5 

Oct. 7-Mar. 22 

n/a 

20' 

lOOj 

n/a 

66 

et al. [1997] 



n/a 

19' 

120J 

n/a 

41 




n/a 

18' 

1701 

n/a 

15 




n/a 

16' 

2501 

n/a 

11 

B. Models 








Schoeberl 

1988-9 

Jan. 3-Feb. 10 

450a 

n/a 

110b 

0.31±0.12a 

n/a 

etal. [1992] 



360a 

n/a 

240b 

0.24±0.06a 

n/a 

Rosenfield 

1988-9 

Nov. 1-Jan. 22k 

500U 

20.2*4 

50b 

1.20c 

34c 

etal [1994] 

& 199 1-2 Nov. 1 -Jan. 22k 

465*4 

n/a 

100b 

1.02c 

n/a 



Nov. 1-Jan. 22 k 

430U 

17.1*4 

140- 150b 

O 

oo 

o 

35c 



Nov. 1-Jan. 22 k 

390'. 1 

1 5.4*4 

200-2 10b 0.72c 

31c 



Jan. 22k-Mar. 21 

445U 

15.5*4 

50b 

0.95c 

81c 



Jan. 22k-Mar. 21 

420'. 1 

n/a 

100b 

0.78c 

n/a 



Jan. 22k-Mar. 21 

395'.i 

1 3.6*>i 

140- 150b 0.60c 

60c 



Jan. 22 k -Mar. 21 

355'.i 

12.2*4 

200-2 10b o.60c 

55c 

Strahan 

n/a 

Oct. 

460 m 

n/a 

250" 

1.2m 

n.a 

etal. [1994] 


Dec. 

415m 

n/a 

250" 

0.75m 

n/a 



Dec. 

470m 

n/a 

1 5 0 n 

0.9 m 

n/a 
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Jan. 

450 m 

n/a 

150" 

0.8m 

n/a 

Feb. 

435m 

n/a 

150" 

0.6m 

n/a 

Mar. 

435m 

n/a 

150" 

0.2m 

n/a 


Lucic 1991-2 Dec. 7-Mar. 15 

475 

n/a 

70" 

0.65-0.9 

n/a 

etal. [1999] 

400 

n/a 

250" 

0.55° 

n/a 


“Estimated from Schoeberl et al. [1992], Figure 4a-parts 3 & 4. 

b Estimated from N20:9 average profile for Jan. 3-Feb. 10, 1989 in Schoeberl et al. [1992], Figure 
4a-part 1. 

c Calculated from 0 (or z) differences. 

dEstimated from observed N2O:0 profiles in Strahan et al. [1994], Figure 3. 

^Estimated from N20:z profile average for Jan. -Feb. 1989 in Loewenstein et al. [1990], Figure 4 . 
fEstimated from N2O:0 profile in Abrams et al. [1996a], Figure 2. 
gEstimated from N20:CH4 correlation fit (this paper). 
bEstimated from difference in 0 (or z) divided by 150 days. 

'0 (or z) of ending profile. 

jEstimated from N20:CCl2F2 correlation in Hartmann et al. [1997], Figure 5. 

January 22 corresponds to middle date of flights reported by Schoeberl et al. [1992]. 

■Estimated from Rosenfield et al. [1994], Figs. 4-5. 

^Estimated from Strahan et al. [1994], Figure 5. 

"Estimated from SKYHI N2O:0 profiles in Strahan et al. [1994], Figure 2. 

°Estimated from Lucic et al. [1999], Figure 6. 


0 — binned data 
binned (inverted) ddtd 









N 2 0/ppbv CH 4 / ppbv 








o 

in 


o 


3 


100 150 200 250 300 0 50 100 150 200 250 300 

N 2 0/ ppbv N 2 0/ppbv 











N 2 0/ ppbv 







0 50 100 -50 0 

Doy of yeor 2000 Day of year 2000 






800 



N 2 0/ppbv CH 4 /ppbv 










Xop J94S 


xop j»d x/iu»3««a 












50 ppbv 
100 ppbv 
150 ppbv 
200 ppbv 
250 ppbv 



Day of year 2000 






